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ABSTRACT 



The Mass loss of Evolved StarS (MESS) sample offers a selection of 78 asymptotic giant branch (AGB) stars and red supergiants 
(RSGs) observed with the PACS photometer on-board Herschel at 70 /um and 160//m. For most of these objects, the dusty AGB 
wind is not spherically symmetric and the wind shape can be subdivided into four classes. In the present paper we concentrate on the 
influence of a companion on the morphology of the stellar wind. Literature was searched to find binaries in the MESS sample, which 
are subsequently linked to their wind-morphology class to assert that the binaries are not distributed equally among the classes. 
In the second part of the paper we concentrate on the circumstellar environment of the two prominent objects R Aqr and W Aql. 
Each shows a characteristic signature of a companion interaction with the stellar wind. For the symbiotic star R Aqr, PACS revealed 
two perfectly opposing arms that in part reflect the previously observed ring-shaped nebula in the optical. However, from the far-IR 
there is evidence that the emitting region is elliptical rather than circular. The outline of the wind of W Aql seems to follow a large 
Archimedean spiral formed by the orbit of the companion but also shows strong indications of an interaction with the interstellar 
medium. We investigated the nature of the companion of W Aql and found that the magnitude of the orbital period supports the size 
of the spiral outline. 

Key words, stars: AGB and post-AGB - binaries: general - stars: winds, outflows - infrared: stars - circumstellar matter 



C\| 1- Introduction velope s differ from spherical sy mmetry while Qlofsson et al.l 

i^i. (1 1988b with CO observations and Kerschba um et al.1 (l2010h in 

L . When low- to intermediate-mass stars evolve off the mam se- ^ far _ TR found ^ ^ ^ detached from ^ star Jn addi _ 

quence and finally up the asymptotic giant branch (AGB) m the tion , asvmmetric s he n s were observed at many wavele ngths: UV 



Hertzsprung-Russell diagram they successively strip off their ( ' ^ ^ 20Q7; Sahai & Chronopoulos 2010), 21 cm Hi 

envelope and lose most of their mass in the form of a dust- ( Gardan et al 20Q6; Matthews et al 2008? 2QU) co ( 



>< 

C3 enriched stellar wind. These dusty outflows are eventually dissi- iK^neXlud^ 

pated into the interstellar medium (ISM) where they contribute ( Ueta et aJ 2QQ6 200g; Gdse et al 20l0; Lad]al et aJ 2Ql0; 

essentially to the chemical enrichment of our Galaxy since the Jorissen et al 2QU Cqx et al 2Ql2) 

vast majority of stars evolve through this stage. 0ne major ^ for a etries is the nce of a 

In the simp est manner the dusty winds expand with a ye- co nion t0 the mass . losing AG B star. The impact of bi- 

locity of 5-l5kms- as spherically symmetric wind bubbles narf Qn ^ sh . of ^ AGB mass 1qss d ds Qn ^ 

that push away the surrounding ISM and lead to e xtended or- h{ separation. In close b inaries, where the AGB star 

cumstellar envelopes (CSEs) around AGB stars ^Young et al.| fiUs ks Roche lobe J Paczvnski | P197D. mass transfer is usu- 

II99J. However, observ ations of post - AGB objects in t faein - ally considered t0 be dynamically un stable (e.g., & \lM 

frared (IR) and optical (|Me,xneret alJ Jl^ Peta et al. |M) \ s„ h Z2 et a lJ ll997t fesen 2004), provided that the AGB 

showed that at least at the end of the AGB phase the en- staf Js more massiye than itg companion; which results in the 

* Herschel is an ES A space observatory with science instruments pro- formation of a common envelope when the accreting star can- 

vided by European-led Principal Investigator consortia and with impor- not adjust fast enough to accommodate the arriving matter (e.g . 



tant participation from NASA. " iLivio & Soken Il988t iTaam & B odenheiiiien ll989USokeJ 1 199 it 
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Hjellming & Taam 1199 It iRasio & Livid fl 996). This common 
envelope is eventually ejected and the system then looks very 
much like a planetary nebula (most likely bipolar) with a binary 
(or even coalesced-binary) nucleus (e.g.. iMiszalski et aDl2012l 
and references therein). 

A less violent variant of this scenario involves a wider system 
where it is not the stellar photosphe re, but rather a dense AGB 
wind that is fi lling the Roche lobe. iMohamed & Podsi adlowski 
(120071 1201 lb dubbed this situation "wind Roche-lobe over- 
flow". The mass transfer is non conservative and matter es- 
capes from the system through an Archimedean spiral, a sit- 
uation already encountered in hydrodyn amical simulations o f 
mass transfer in de t ached syste ms by e.g. ISawada et all dl986). 
The uns & Jorissenl (119931) and iMastrodemos & Morrisl dl998. 
19991). where this feature was called "s pira l arm" (see also 
Huggins et alj 120091: iKim & Taarr] l2012aB l3: IWang & Willsonl 
20121) . Observationally, such lar ge-scale spirals were ob- 



served around LL Peg (AFGL 3068; Mauron & Huggins 2 006 
Mor ris et al.l 120061) . RW L Mi (CIT 6: iDinh-V.-Trung & Liml 
2009b, and recently R Scl dMaercker et al.l 120 12b. and als o 
in a bruised version around o Cet (Mira; iMayer et al.l 1201 ll) . 
Estimates of the orbital separation are available for three 
(AFGL 3068, CIT 6, and Mira) of these four systems, flagging 
them as wide bi naries given their projected separations of 109 , 
70, and 55 AU dMorris et ai1l2006t IDinh-V.-Trung & Limll2009t 
iKarovska et al.ll 19971) . In the framework of this paper, which in- 
terpret s Herschel observa tions, it is important to report the find- 
ing bv lEdgaret all (12008). The authors suggested that the spiral 
shock associated with mass-losing AGB stars in binary systems 
may anneal the amorphous grains to form crystals, as is ofte n 
seen in the discs around post- AGB binaries (iGielen et al.ll201ll) . 

Finally, one should mention the freq uent occurren ce of bipo- 
lar jet s in binary AGB sta rs (e.g., Huggins 2007), such as 
o Cet dMeaburn et al.l 120091) . R Aqr jWallerstein & Greensteinl 
1980; Kafato s et al.1 fl989). tt 1 Gru (|Sahai| ll992|), and V Hya 
( Hirano et aDl2004 



In this paper, we examine a sub- sample of the 78 AGB stars 
from the MESS (Mass loss of Evolved StarS; IGroenewegen et al.l 
l201ll) guaranteed-time key program (GTKP) for Herschel and 
give new observational insights into this topic. In Sect. 13.21 we 
give a short statistical overview of the morphology of circumstel- 
lar envelopes of binaries in the sample. The focus is then set on 
two prominent examples: R Aqr and W Aql. Each of them shows 
remarkable structures in the CSE that point to various forms of 
interaction. 



2. Observations and the MESS sample 

2.1. Observations 

The observations presented here were carried out with 
the Photodetector A rray Camera and Spectrometer (PACS; 
iPoglitsch et aDl2Q10h on-board Herschel in 2009-2011. PACS 
is operating in the far-IR wavelength range in three channels, 
70yum (blue channel), 100 yum (green channel), and 160 yum (red 
channel) two of which are used simultaneously. Usually, these 
are the blue and red bands. The pixel size in the two channels is 
3'/2 x 3'/2 in the blue band and 6'.'4 x 6'.'4 in the red band with 
a typical field of view of 5' x 5', but the images presented here 
are oversampled by a factor of 3.2, resulting in a resolution of 
1" pix" 1 in the blue and 2" pix" 1 in the red band. This image 
resolution has to be compared to the FWHM of 5 '.'7 and 1 1'/4 at 
these wavelengths. 



The data processing was performed using the standard re- 
duction steps in the Herschel interactive processing environ- 
ment (HIPE) and the map-makin g program Scanamo rphos. Both 
are described in detail by iGroenewegen et al.l (1201 ll) and Rous- 
sel (in prep.). In addition, deconvolution was applied to all 
images. A detailed a nalysis of this process can be found in 
lOttensamer et al.l (1201 ll) . 

2.2. MESS program 



The MESS GTKP I Groenewegen et al.1 1201 ll) performed with 
the Herschel Space Observatory dPilbratt et al 11201 Ol) strives to 
investigate the CSE and the mass loss history of evolved stel- 
lar objects, benefits from the high spatia l resolution of Her- 
schel's PACS and SPIRE instruments (IPoglitsch et al.l I201QI: 
iGriffin et al.ll2010l) . 

Earlier results of this program included the detections of 
multiple shells to gether with a bow shock around CW Leo 
(IRC +10216) bv ILadial et all (l2010h and IDecin etaD (l201ll) 
and of detached shells a ro und AQ And, U An t and TT Cyg 
dKerschbaum et al. 2010). iJorissen et al.l (l201ll) examined the 
envelopes of TX Psc and X Her. The envelopes are massively 
shaped by the interaction of the stellar winds with the oncoming 
IS M due to their high s pace velocity. 

Ma yer et al.l (l201ll) discussed the circumstellar environment 
of o Cet and found it to be shaped by a combination of several 
processes: wind-ISM interaction, a bipolar outflow, and in the 
interior, the orbiting motion of the companion that forms a spiral 
that appears to be broken. The broken appearance results from 
the bipolar outflow that pierces the dust shells. 

2.3. MESS sample categories 

For this study on binary-induced structures, we concentrate on 
the 78 AGB stars and RSGs in the MESS sample. The compo- 
sition of this sample is well-balanced since it contains all major 
spectral classes of late-type stars. More specifically, the sample 
consists of 37 C-rich, 32 O-rich stars and 9 S- stars. 

Recently, ICox et al.1 ([2012) reviewed the properties of the 
whole MESS sample with the aim to classify the different CSE 
morphologies and to relate them to the influence of the interact- 
ing ISM. In short, five main morphology classes were identified: 

- the fermata type, which is characterised by a strong bow 
shock in the direction of the space motion of the star; 

- the eye type, or double-fermata, which has two elliptically 
shaped density-enhanced regions; 

- the ring type, which is in general a detached shell (see e.g. 
iKerschbaum et al.ll20Tol) : 

- the irregular type, which has a completely asymmetric mor- 
phology. This can be due, for example, to very strong insta- 
bilities that destroy the bow shock or to strong interactions 
with a companion; 

- a point source, which shows no resolved CSE at all. 



A detailed discussion is provided by Co x et al ] (120121) and will 
not be repeated here. The strong indications of ISM interaction 
and the wide variety of morphologies motivated us to investi- 
gate the possible role of companions that shape the CSE in more 
detail. A discussion of that role based on the PACS images is 
the topic of the present paper. We refer to the above-mentioned 
categories throughout. 
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Table 1. Properties of the suspected binaries in the MESS sample. The table is subdivided into two pa r ts. Th e upper part lists the physical binaries, the lower part the visual binaries without 
indicating physical association. The sp ace motion v LS r is c alculated following John son & Soderbloml (Q987). The adopted solar motion to convert heliocentric motion into LSR motion is 
(U,V,W)o = (11.10, 12.24, 7.25)kms _1 dSchonrich e t al. 2010). The objects are ordered by increasing distance. 



Object 


Morphology 


Spec. Type 


d 


M 

riO -7 Mr, vr _1 l 


v w 

[kms -1 ] 


''LSR 

[kms -1 ] 


Binary detection 


References 


o Cet 


Fermata & Irregular M7 Hie 


91 + 10 


2.5 


8.1 


107.7 


1 t~\l 1 f T~» A AAO 

visual: p = 076, PA « 90 


(1),(2),(3) 


Aps 


Fermata 


M6.5 III 


113 + 6 


1.1 


4.5 


34.2 


Hipparcos variability-induced mover (VIM) 


(1),(2),(4),(5) 


EP Aqr 


Fermata 


M8 IIIv 


114 + 8 


3.1 


11.5 


40.2 


complex CO profile 


(l),(2),(6),(7) 


1 

n Gru 


Irregular 


S6e 


163 + 20 


8.5 


30.0 


11.9 


jets, visual: p = 2745, PA = 200.4 


(l),(2),(8) 


1 r\ • 

o l On 


Irregular 


M3 


200 + 28 


< 0.4 




38.8 


UV excess 


(9),(l0),(ll) 


R Aqr 


Irregular 


M7 Illpev 


214 


0.6 


16.7 


39.7 


symbiotic 


(I2),(l3),(l4),(l5),(l6),(l7) 


Y Lyn 


Point source 


M6 SIb-II 


253 + 61 


5.6 


5.4 


13.9 


Hipparcos VIM 


(l),(9),(5) 


TW Hor 


Point source 


C 


322 + 38 


0.24 


7.5 


52.4 


UV excess 


(I),(l8),(l9),(l6) 


W Aql 


Fermata 


S6.6 


340 


130 


20.0 


21.6 


composite spectrum/CO profile, 
visual: p = 0747, PA = 209° 


(20),(2),(1),(21),(18),(22) 


R Scl 


Fermata & Rin 


g c 


370 


16 


17.0 


66.3 


spiral arm 


(2),(l)(23) 


VY UMa 


Eye 


C 


380 + 51 


0.7 


7.9 


28.2 


T TV pYppss 


n^i (4-^ ad} no^i 

V 1 ) A v A^^J A 1 u / 


TX Cam 


Point source 


M8.5 


380 


65 


21.2 




spiral arm 


(2),(19) 


UCam 


Eye 


C 


430 


10 


20.6 


9.8 


Hipparcos DMSA/C: p = 0717, PA = 79° 


(25),(26),(27),(1) 


VEri 


Point source 


M5/M6IV 


439 + 133 


1.5 


13 


39.8 


UV excess 


(1),(26),(24),(10) 


RWLMi 


Point source 


C 


440 


59 


20.8 




spiral arm 


(6),(2),(28) 


VY CMa 


Point source 


M3/M4 III 


562 


2800 


46.5 


42.9 


UV excess 


(26),(2),(1),(10) 


VHya 


Point source 


C 


600 


610 


30 


27.9 


rapid rotation 


(14),(2),(1),(29) 


LL Peg 


Point source 




980 


310 


16.0 




spiral arm, visual: p = 0711 


(2),(30),(31) 


R Dor 


Point source 


M8 Hie 


55 + 3 


6.1 


6.0 


39.4 


visual: p = 3273 


(32),(33),(34) 


W Hya 


Ring 


C 


104 + 12 


0.78 


7.5 


48.2 


visual: p = 7475 


(1),(2),(34) 


R Hya 


Fermata 


M7 Hie 


124 +11 


1.6 


12.5 


25.3 


visual: p = 2172 


(I) (2) (34) 


RCas 


Fermata 


M7 Hie 


127 + 16 


4.0 


13.5 


44.0 


visual: p = 2778 


(D,(2),(34) 


RCrt 


Eye 


M7 III 


261 + 65 


6.9 


10.8 


23.8 


visual: p = 6574 


(D,(2),(34) 


VX Sgr 


Point source 


M5/M6 III 


262+ 187 


610 


24.3 


8.5 


visual: p = 077 


(1),(2),(34) 


UU Aur 


Fermata 


C 


341 


2.7 


11.0 


17.9 


visual: p = 1872 


(32),(2),(34) 


VPav 


Eye 


C 


370 + 73 


3.4 


16.0 


24.0 


visual: p = 1570 


(1),(14),(24),(34) 


S Set 


Ring 


C 


386 + 85 


1.4 


8.0 


17.5 


visual: p = 1473 


(D,(2),(34) 


TT Cyg 


Ring 


C 


436 


50 


12.6 


38.3 


visual: p = 7170 


(32),(35),(36),(34) 



Refer ences. (D:lvan Leeuwenl d2007h. r2):lDe Beck etan d2010h. G ): lKarovska eTall dl997h. (4 ): ISchoier & Olofssonl j200lh. ffl : IPourbaix et all d2003h r 6):lKnapp et alJdl998h. r7):IWinters et all 
J2007L (8): ISahail d!992l), (9 ): iGroeneweg en & de Jond d!998l) (10): ISahai et all d2008l) (11): lAke & Johnsonl d!988h. (1 2): iKamohara et al I d2010h . (13): iMatthew s & Reid j 2007 [) . 
(14): [Groenewegen et all d2002l). (15): IWillson et all d!981l) . (1 6): iHinkle et all d!989l) . (17): iGromadzki & Mikolai ewskal d2009h. (18): ICulver et al.1 ( 119741). (19): ICastro-Carrizo et all d2010h 
(20): [ Guandalin i & Bussol d2008l). (21): iHerbid d!965h. (22): iRamstedt et all 1201 In. (23): iMaercker et all d2012h. (24): lLoupetal.1 d!993h . (25): iKnapp et all d2003l) . (26): I Young et al.1 
dl99 3[). (27): iLindavist et all d!999h. (28): iDinh-V.-Trung & Liml (120091). (29): Barnbau m et al I d!995h . (30): iMauron & Hugginsl d2006l) . (31): iMorris et al.1 d2006l) . (32): ICox et all d2012h . 
(33): lQlofsson et alJfc002h T734): lProust et al.ldl981h . (35): l01ofssonl ([T996). (36): blofsson et al.ld2000h 

Notes. Although VX Sgr is associated with entry ADS 11078 in lProust et al.l d!981h . the Washington Double Star Catalog ( [http : //ad . usno . navy . mil/wds/wdstext . htmT) does not list any 
measurement for that system, which is thus most likely spurious. 



3. Physical binaries in the MESS sample 

As described above, the FWHM of the PSF is 5'/7 and 11 '.'4 
at 70yum and 160/zm, respectively. However, these resolutions 
are insufficient to resolve the individual components of a binary 
system even for relatively nearby objects. For this reason, the lit- 
erature was thoroughly searched to find other possible evidence 
for binarity among the 78 AGB stars and RSGs in the MESS 
sample. The result of this search is listed in Table [T] resulting in 
18 objects that show strong indications for physical binarity in 
their spectra, proper motion history, or observed morphology. In 
addition, the table also lists the objects of t he MESS sample tha t 
are part of the visual binary catalogue by [Proust et all (Il98ll) . 
but we stress that most of these visual binaries are likely to be 
optical pairs given their large angular separations. 

Table [T] also provides the general and kinematic data of the 
targets. In the penultimate column, the detection method of the 
binarity is listed. 



3.1. Difficulties and limitations of binary detection 

Flagging an AGB star as a binary is a difficult endeavour, mainly 
because the main binary-detection method - variable radial ve- 
locities - is not very efficient in the case of AGB stars (no con- 
vincing case of a spectroscopic binary involving a genuine AGB 
staiQ is currently known). This lack of efficiency results from the 
fact that orbital radial- velocity variations are confused by the in- 
trinsic variations a s sociated with the atmosphe ric dynamics (e.g., 
iHinkle et al.lll982ll2QQ2l:lNowotnv et al]l2010h . especially since 
the large radii of AGB stars impose their orbits to be rather wide, 
and hence the velocity amplitude to be small. 

Other methods of binary detection should thus be used in the 
case of AGB stars, of which we list: 

- detection of the companion itself (visually, interferometri- 
cally or through lunar occultations); 

- composite nature of the spectrum and/or symbiotic activity; 

- X-ray emission; 

- shallow light curve of a Mira variable, with wide and flat 
minima; 

- light variations associated with orbital motion (eclipses, Al- 
gol); 

- asymmetries in circumstellar envelopes (especially spiral 
features, or double CO line profile); 

- proper-motion variations. 

A detailed description of these methods ( with a ppro- 
priate references) may be found in Jorissen I (E004) and 
iJorissen & Frankowskil (I2008I) . The selection of putative bina- 
ries among the MESS targets is based on these different meth- 
ods, as indicated in the penultimate column of Table [T] The con- 
fidence level for a target to be a possible binary varies from one 
target to another. In the case of visual binaries with a small 
angular separation (o Cet, n l Gru, W Aql), there can be no 
doubt; there is very little doubt in the case of UV-excess de- 
tection (o l Ori, VY UMa, TW Hor, V Eri), but the confidence 
is lowest for the double CO profile (EP Aqr). To account for 
these profiles, a bipolar geometry accompanying a Keplerian 
circumbinary disc has been proposed for instance for BM Gem 
(Kahane & Jura 1996), whose companion was later detected by 
Izumiur a et al.l (120081) . But a double wind (fast and slow), as 

1 a star is defined as an AGB star if it is a Mira variable, carbon- rich, 
or exhibits spectral lines from the unstable element technetium 



Table 2. Distribution of morphologies in the MESS sample. 



Type 




Number 




Number 


Physical binaries 








within 500 pc 


within 500 pc 


Fermata 


25 


[32.1%] 


23 


[39.0%] 


5/23 


[21.7%] 


Eye 


7 


[9.0 %] 


7 


[11.9%] 


2/7 


[28.6 %] 


Ring 


15 


[19.2%] 


13 


[22.0 %] 


1/13 


[7.7 %] 


Irregular 


8 


[10.3 %] 


7 


[11.9%] 


4/7 


[57.1 %] 


Point source 


30 


[38.5 %] 


15 


[25.4 %] 


5/15 


[33.3 %] 


Total 


78 




59 




15/59 


[25.4 %] 



Notes. Seven objects, of which two are binaries (o Cet and R Scl), show 
multiple morphologies and are counted for t wo types. Of the se seven 
objects, six are within 500 pc. For details see lCox et al .1 (12012). 



originally suggested by Knapp et al. (1998) for instance for the 
MESS target EP Aqr is an alternative to the binary model. The 
recent CO map of the circumstellar environment of that object by 
Nakashima ( 2006) casts doubt on the binary nature of that star, 
since it does not reveal any signature of a circumbinary disc. In 
those cases, other binary diagnostics, such as the effect of bina- 
rity on the shaping of the CSE can dispel doubts. 

3.2. Distribution of binaries according to morphological 
categories 

In this section, we investigate the statistical frequency of bina- 
ries among the different morphological categories to see whether 
some categories are more often associated with binaries than oth- 
ers. 

Before embarking on this comparison, we emphasize, however, 
that our ability to detect complex structures around AGB stars 
depends upon their distance, since it is obviously easier for 
near by objects. Thi s is reflected in the statistics derived from 
the lCox et al.l (120121) morphologies: the number of point sources 
drops from 30/78 (38.5 % of the total sample) to 15/59 (25.4 %) 
for the d<500pc objects, as seen in Table O For this reason, 
the binary frequency among the different morphologies is only 
investigated for the 59 objects within 500 pc from the Sun. 

Only three among the 18 binaries are farther away than 
500 pc, which yields a binary frequency of 25.4 % in the MESS 
d<500pc sample. T his value is not so d ifferent from the 30 % of 
binaries detected by Ivan Winckell (120031) among post- AGB stars. 

Furthermore, seven objects of the MESS sample show mul- 
tiple morphologies, and two of them are binary stars: o Cet, 
whose extended s tructure was already discussed in detail by 
iMaver et al] (1201 ll) . and R Scl dMaercker et al. 11201 2b . 

The most common morphology type within 500 pc of the 
sample is the fermata: 23 out of the 59 objects (39.0%) show 
this kind of strong AGB wind - ISM interaction. Of these 23 
objects, 5 (or 21.7%) are physical binaries. Do these five bi- 
nary stars have a fermata morphology, which would suggest that 
the main shaping agent is indeed the interaction of the fast mov- 
ing AGB wind with the ISM, which would mean in turn that 
the binary motion is negligible? In fact, the classification is not 
as clear-cut as it might seem. For instance, the binary nature 
of EP Aqr is subject to debate, as discussed in Sect. 13.11 and 
should thus not be included in the present discussion. For o Cet 
and R Scl , ferma ta is not the sole morphological type, because 
ICox et al.1 (120121) also associate them with types "irregular" and 
"ring", respectively. Since o Cet is the nearest binary target, it 
is not really surprising that this system shows the most conspic- 
uous impact of the binary motion, in the form of broken arcs 
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and of a symmetry axis tilted with respect to the space motion 
(iMayer et al .1201 II) . 6 Aps and W Aql, like o Cet, are other cases 
where signatures of binary motion and ISM interaction seem to 
co-exist in the CSE appearance: we show in Sect. 14.21 that an 
Archimedean spiral fits the CSE of W Aql well, even though the 
spiral also shows some signs of interaction with the ISM. In a 
forthcoming paper, we will demonstrate that the CSE of 6 Aps 
shows indications for a jet bumping into the wind — ISM bow 
shock. Thus, there is not one single pure fermata morphology 
among the binary stars. 

The eye type is a comparatively small group of seven objects 
(representing 11.9 % of the d<500pc sample) two of which are 
binaries (28.6%). Given the small-number statistics, no conclu- 
sion can be drawn for this category. 

The ring category contains 13 members within 500 pc 
(22.0%). Ring objects in general are slowly moving through 
the ISM, and most of them are carbon stars. In some cases, how- 
ever, a bow shock forms beyond the ring structure, supporting 
the assumption that the rings are rather due to wind- wind inter- 
action, probably asso ciated to the tra nsition from an O-rich to 
a C-rich atmosphere (iCox et al. I l2012h . If for one reason or an- 
other, the wind from the C-star becomes faster than the wind 
formerly present around the O-star (for instance because of the 
change in opacity of the carbon grains with respect to oxygen- 
based grains, or because of the increase in luminosity accom- 
panying the thermal pulse responsible for the carbon synthesis), 
the fast wind will cat ch the old slower wind and interact with it 
(iMattsson et alj2007l) . probably forming the observed ring. Cur- 
rently, only one (R Scl) of the 13 ring objects sh ows indications 
of bin arity, in the form of an Archimedean spiral (Maer cker et al.1 
120121) . 

The shaping influence of binarity is best revealed by the ir- 
regular type stars. Their circumstellar environment exhibits no 
symmetry and cannot be related to ISM or to wind- wind interac- 
tion. There are seven objects belonging to the irregular category, 
four of which are binaries (57.1 %), which is by far the high- 
est percentage. The number of objects is too small though to 
conclude that binaries are the only cause for an irregular mor- 
phology. The four binary objects in the irregular category are 
o Cet, o l Qri, n l G ru, and R Aqr. o Cet was already discussed by 
Ma yer et all (120111) and for o l Ori there is doubt whether the very 
faint emission seen at the sky noise level with PACS is indeed 
real and associated to the central star (see Fig. 4. in ICox et all 
l20l2l) . The CSE of n l Gru shows an arm structure, which is 
strong evidence for the influence of a companion. This target 
will be discussed together with Aps in a forthcoming paper. 
The IR emission of R Aqr, dominated by two opposing arms, is 
discussed in Sect. 14.11 The remaining stars in the irregular cate- 
gory with n o evidence so far for binarity are AFGL 278, ST Her, 
and V Cyg (ICox et al .1120121) . Despite the small number of ob- 
jects, it is possible to compare the binary frequency in the whole 
sample (pb = 15/59 = 0.254) with that in the "irregular" cat- 
egory (/?B,irr = 4/7 = 0.571) and check the null hypothesis Ho 
that the two frequencies are indeed identical and their difference 
only due to sampling fluctuations. If the null hypothesis is true, 
the expected number of binaries among irregular morphologies 
is 

x = N BM =p B N kr = lJS, (1) 

as may be estimated from the fraction of binaries in the total 
sample (denoted p B ) applied to the number of stars in the irreg- 
ular category (Nfa). The corresponding standard deviation on x 
(denoted cr x ) is computed from its expression for a hypergeo- 
metric distribution, which characterises the variable x, given the 



total number of stars in the sample, N^i, the total number of bina- 
ries, p B N a \\, and the number of stars with irregular morphologies, 

o- x = [N irr p B (l - pb)(Nm - AW/Wdi " 1)] 1/2 = 1.08. (2) 

The reduced, squared difference c 2 between the estimated num- 
ber of binaries with irregular morphologies (x = 1.78) and the 
observed number (x = 4) is 

c 2 = [(x - x)/cr x ] 2 = 4.225, (3) 

and its significance is computed from the x 2 distribution with 
one degree of freedom. The corresponding significance level is 
then 

a = 0.5 P(x 2 > c 2 ) = 0.02. (4) 

Thus, the probability that the difference between the expected 
number of binaries with irregular morphologies (1.78) and the 
observed number (4) is due to random fluctuations amounts to 
only 2 %, and may be considered as insignificant. 

The last type of objects are the point sources, representing 15 
of the 59 objects within 500 pc (25.4 %). Five of them are bina- 
ries (33.3 %), which is about the average percentage of binaries 
in the whole sample (23.7 %). 

Finally, we s tress that the sample selection, as sketched by 
ICox et all(l2012l) . should not interfere with the above statistics, 
since that statistics is based on the distribution of binaries 
among morphological categories defined a posteriori. 

According to these results, a companion to an AGB star 
would most likely produce circumstellar structures that are 
ascribed to the irregular morphology type. However, a closer 
look at the images, which we report in Sect. [4] for R Aqr and 
W Aql, leads to the conclusion that the irregular category hides 
morphologies specific to binary sys tems, such as jets (R Aqr) 
or spiral arms (W Aql and o Cet; IMayer et aD l201l[ for the 
latter). But apparently nothing prevents this binary-induced 
morphology from coexisting with a bow shock that might result 
from interaction with the ISM, which is observed on a larger 
scale, as in the latter two systems. 

4. Case studies of binary AGB stellar winds: R Aqr 
and W Aql 

In this section, the new Herschel-PACS far-IR observations of 
the large-scale environment of R Aqr and W Aql are discussed 
with emphasis on how the binary interplay shapes the AGB wind 
morphology. These objects were selected because of the com- 
plex morphology of their envelope, which suggests a prominent 
role of binary interaction. This section attempts to elucidate that 
role. 

4.1. RAqr 
4.1.1. Background 

R Aqr is one of the best- s tudied and brightest symbiotic sys- 
tems dWhitelock et al .11 983). probably observed for almost 1000 
years, as suggeste d by a nova reco rd in Korean history books 
(1073-1074 A.D.: lYang etaDl2005h . The system consists of an 
M7III Mira variable primar y star with a pulsation period of 3 87 
davs dMattei & Allenll 19791) and a hot companion (lMerrilllll92Tl 
Il950h that is most probably a magnetic white dwarf with M « 
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Fig. 1. Deconvolved PACS images of R Aqr at 70 jim (left) and at 160//m (right). The arrows indicate the direction of space motion and the 
position of the star in 2000 yr. The proper motion values are taken from Kamohara et al. ( 2010). 



0.6 - 1 M Q dNichols et al.1 120071: [Gromadzki & Mikolaiewskal 
2009). The companion was resolved by lHollis et al.l (1 19971) with 
a separation of 55 ± 2 mas from the primary (1 1.8 AU at 214 pc) 
at a position angle of 18° ± 2°. 

The orbi ta l per iod of 44 yr was first determined by 
Will son et al.l (Il98ll) from eclipses in the light curves in 
1928-1934 and 1974-1980 and w as refined later to 43.6 yr 
( Gromadzki & Mikolajewska 2009). The orientation of the orbit 
during the eclipse is uncertain, h owever. From the duration of 
the eclipse, Makinen et al. ( 2004) claimed that it occurs at apas- 
tron, while iMcIntosh & Rustanl (120071) suggested that it happens 
during periastron passage. The latter authors also obtained an 
orbital period of 34.6 yr, significantly shorter than the other esti- 
mates. 

The orbit is highly eccentric, with e « 0.6-0.8, and i t is see n 
nearly edge-on wi t h i = 70° as noted by lSolf & Ulrichl (TT985). 
IWallersteinl([T986l) . lHinkle et al.l(ll989l) . and lHollis et aD(ll997h . 

The distance to the system is fairly well estab lished and 
ranges from 363 p c (Hipparcos; Ivan Leeuwenl l2QQ7h to 214 pc 
(SiO maser spots; iKamohara et al.l 120101) . wit h an intermedi- 
ate v alue at 240 pc (period-luminosity relation; Iwhitelock et al.l 
2008]). T hroughout this d i scussi on, the distance of 214 pc is 
used since IKamohara et al.l (120101) were also able to derive new 
proper motion values with a small error bar (ji a cos S = 32.2 ± 
0.8 mas yr -1 and^ = -29.5 ± 0.7 mas yr" 1 ). 

In the optical, the system is surrounded by two jets and 
a nebular structure with the shape of a double hourglass 
(ISolf & Ulrichlll985h. The NE optical jet (PA = 46°) was dis- 
covered first (Wa llerstein & Greensteinl H980I) at a distance of 
about 6" from the binary system, while the sec ond (SW) jet 
was discovered as knots in the radio emission (Kafatoset aD 
198 9 ). Subsequently, the jets were intensively studied (e.g . 



Holl is et al.lll99QL 119911: lyTotti et alJll987b iMakinen et al1l2004l: 
Nichols et al.ll2007l) . iKellogg et al. (200% observed both jets in 
X-ray and radio emission and found the NE jet to be curved and 
dominated by a bright spot 8 " from the star. The br ight spot is 
tilted away from the observer dNichols & Slavinll2009l) . The SW 
jet, with PA = 211°, is not directl y opposite, with a b right spot 
26" away from t he binary system (Kellog g et al.ll2007l) . 

Interestingly, Makine net al.l d2004) found that parts of the 
NE jet broke apart within a timespan of 7 yr probably due to a 
collision with a more massive cloud. 

In addition to the jet, the circumst ellar environment fea tures 
an hourglass structure (discovered by ISolf & Ulrichl IT9 85). ex- 
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Fig. 2. SED of R Aqr described by single-band photometry (filled 
squares) compiled from the literature as well as the available ISO-SWS 
(blue line) and -LWS (cyan line) spectra. Overplotted are the synthetic 
spectra based on models for the photosphere (MARCS; g reen lin e) and 
the combined photosphere + wind (red line; see Section I4.1.2I for de- 
tails). In addition, aperture photometry of the PACS 70 iim and 160//m 
bands (triangles) is included. There are three data points for each of 
these measurements. From bottom to top: aperture photometry (i) of the 
eastern arm, (ii) of the central source, and (iii) of the whole IR-emitting 
region. The dotted line represents the model for the dust emission of 
the eastern arm. 



tending to about 2'. It is centred on the binary system and might 
be linked to a possible previous explosion on the white dwarf 
660 yr ago. The structure of the nebula consists of a ring with 
a ra dius of 42" inclined by 18° with respect to the plane of the 
sky. iHollis et all (Il989b compared 6-cm radio continuum emis- 
sion with observations in the optical and found them to be mor- 
phologically similar and spatially related. 



4.1 .2. Characterising R Aqr and its inner wind region 

To understand Herschel's observations of the outer regions of the 
environment of R Aqr, it is important to derive characteristics of 
the star, its present-day mass loss and the region where the out- 
flow is triggered. For this reason, the literature was searched 
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Table 3. Proper motion values of R Aqr and W Aql from various catalogues sorted according to their mean epoch. \i a * stands for \i a cos 6. 



R Aqr 


GC (1) 


SAO (2) 


Tycho-2 (3) 


UCAC3 (4) 


PPMXL (5) 


USNO (6) 


HIP2 (7) 


ASCC-2.5 (8) 


SiO maser (9) 




i on'* i 






1QQ1 07 

YyyY.L 1 


1 QQ/I £(\ 


zUUU.U 








fi a * [mas/yr] 


27 


29 ±5 


33 + 1.4 


30.8 + 1.6 


28.6+1.7 


32 


33.00 + 1.53 


28.95 + 1.52 


32.2 + 0.8 


ld 5 [mas/yr] 


-20 


-21 +5 


-32.6+ 1.2 


-31.5 + 2.1 


-31.8 + 1.6 


-34 


-25.73 + 1.30 


-32.61 + 1.85 


-29.5 + 0.7 


WAql 


GC^ 


SAO< 2 > 


YZC (io) 


W550 (11) 


CPIRSS (12) 


UCAC3 (4) 


CMC (13) 


ASCC-2.5 (7) 


PPMXL (6) 




1913.2 


1913.2 


1933.60 


1968.55 




1981.86 


1993.68 




2002.26 


yi/ a * [mas/yr] 


27 


30 + 20 


26 


30 


33 


-5.2 + 48.1 


15.9 + 1.8 


15.06 + 1.70 


-21.0 + 20.8 


fi s [mas/yr] 


7 


8 + 22 





10 


8 


17.6 + 48.1 


-0.9 + 2.0 


-0.58 + 2.00 


- 36.5 + 20.8 



References. (1): General Catalogue o f Stars dBossI |1937 | ). (2): Smi thsonian Astrophysical Observatory Star Catalog 

dSmithsonian Astrophvsical qbservatorv StafU ll966). (3): Tvcho-2 (lH0g et all [2000h . (4): Third U.S. Naval Observatory CCD Astrograph 
Catal og dZacharias et alj2010L based on Hipparcos and Tycho- 2 for bright and SP M and SuperCosmos for faint stars), (5): PPMXL dRoeser et al.l 
120101 . combining US NO-B1.0 and 2MAS S), (6): USNO-B1.0 dMonet et al.ll2003l taken from 7435 Schmidt p lates over the last 50 yr), (7): New 
Hipparcos Reduction ( van Leeuwen 2007 ), (8): All-Sky Compiled Catalogue-2.5 ( Kharchenko & Roeser 2009, merging e.g. Hipparcos, Tycho-1, 
Tycho-2 and PPM), (9): Kamo hara et"alH2010h. (10): Yale Zone Catalogue (Fallonl l 1983k (11): Washington 550, Results of observations with 
the six-inch transit circle 1963 -1971 (Hu ghes & Scott|[l982h . (12): The U.S. Nava l Observatory Catalog of Positions of In frared Stellar Sources 
dHindslev & Harrington 1991), (13): Carlsberg Meridian Catalogs Number 1-11 (Carlsberg Meridian Catalogues Staff 1999) 



Table 4. Input parameters used for modelling the atmosphere and wind 
of R Aqr. 



luminosity L+ [L ] 
effective temperature [K] 
mass M* [M Q ] 

surface gravity log(g*) [cm s" 2 ] 
metallicity Z [Z©] 
C/O ratio 

mass-loss rate M [M© yr" 1 ] 



4780 
2800 
1 

-0.5 
1 

0.55 
7 x 10" 7 



to derive a spectral energy distribution (SED) of R Aqr. Photo- 
metric data were adopted from the ASCC-2.5 ca talogue (B, V 
bands) , the 2MASS catalogue (J,H,K bands), iFouque et al.l 
(119921) [L, M bands], the Infrared Astronomical Satellite (IRAS; 
12, 25, 60 & 100 /mi), and the AKARI satellite (65, 90, 140 & 
160/mi). Valuable information is provided by the spectra ob- 
tained by the Infrared Space Observatory [ISO; 2-45 /mi (SWS) 
and 43-160yum (LWS)]. Figure [2] shows the resulting SED of 
R Aqr. The short- wavelength band (SWS) of ISO is dominated 
by pronounced emission of an optically thin dust shell that con- 
tains amorphous silicate grains. No signs of crystalline dust fea- 
tures have been found in R Aqr so far. In ISO's long- wavelength 
range (LWS), the spectrum shows no prominent features. The 
ISO spectrograph had a wavelength-dependent field of view of 
« 25" x 25 ", thus not covering the outer regions observed by 
Herschel. To derive the Herschel-PACS fluxes of the central 
star, aperture photometry with a radius of 10" was applied. The 
resulting fluxes are F v jq = 20.980 Jy and F Vj i6o = 2.569 Jy at 
70//m and 160//m, respectively, and coincide well with the pre- 
dictions of the combined photosphere and wind model (to be 
described below). 

Three PACS fluxes are plotted in Fig. |2] The first (highest) 
one corresponds to the overall IR emission (aperture photom- 
etry on a circle of radius 50" with resulting fluxes of F v jo = 
28.220Jy and F v? i 60 = 4.320Jy at 70 //m and 160 /mi, respec- 
tively). The intermediate flux value corresponds to the above 
described measurement of the central source, while the third 
(weakest) one corresponds to the flux of the eastern arm. An 
elliptical aperture with a size of 40" x 20" tilted by 45° was cho- 
sen to fully cover the arm. The resulting fluxes for the arm are 
F Vf70 = 0.481 Jy and F y , 160 = 0.288 Jy at 70 a nd 16 jam, re- 
spectively. As discussed in more detail in Sect. 14.1.31 this flux 



measurement in the CSE of R Aqr may be interpreted as astro- 
nomical silicate dust emission. Adopting absorption coefficients 
2absjo = 3 x 10" 3 and Qahs i6o = 5.9 x 10" 4 at 70 fim and 160 /mi, 
respectively (lDraine| [l985), the 10 jam and 160 jam fluxes ob- 
served at the eastern tip o f the eastern arm th en correspond to 
a temperature of 37 K (see Joriss en et al .1 1201 1L for details). Be- 
cause there is a temperature gradient in the arm (see Sect. 14.1.3b , 
the average temperature of the arm is ^ 32 K. The corresponding 
dust-emission spectrum - computed via vnB v (T - 32K)g a b s (v)/, 
where / is a scaling factor - is also shown in Fig. [2] 

To characterize R Aqr, we modelled the atmosphere and the 
inner wind region using methods of proven reliability. The 
model parameters, listed in Table ffl were chosen according to 
values compiled from the literature where available . With an an- 
gular diameter of 6* = 0V015 ( van Belle et a 1. 1997), a distance 
of D = 214 pc dKamohara et al.ll2010b , and assuming the star to 
have a mass of 1 M©, one can estimate l ogp?*). Together with 
the e ffective temperature adopted from IContini & Fo rmiggini 
(2003), the luminosity of R Aqr can be derived. Metallicity 
and C/O ratio were assum ed to be solar- l ike fo llowing the abun- 
dance pattern derived by ICaffau et al] (I2008L 12009b . For the 
stellar atmosphe re, we compu t ed a C OMARCS model atmo- 
sphere following Aringer et al. (2009). This was then comple- 
mented with a model for the stationary outflow containing amor- 
phous silicate dust parti cles, calculated as described in detail by 
iFerrarotti & Gaill (120061) . including recent improvements of the 
wind code (Gail et al., in prep.). Based on this combined model, 
we carried out detailed radiative transfer calculations sim ilar to 
previous works (I Aringer et al . 2009; N owotny et al.ll20lH) . 

This led to a synthetic spectrum that simulates an unresolved 
object with spectral contributions from the photosphere and the 
inner- wind region (but not from any extended circumstellar ma- 
terial farther than the eastern arm of R Aqr), just as can be as- 
sumed for the observational data in Fig. |2] The resulting model 
spectra were normalised to the flux of the ISO- SWS spectrum 
at 3.8 yum for a direct comparison. Overplotted in the figure is 
the best-fitting model after iterating the mass-loss rate. While 
the COMARCS model is able to reproduce the SED of R Aqr in 
the visual and NIR reasonably well, an additional emission due 
to the circumstellar dust is clearly needed for a proper descrip- 
tion beyond « 4 jim. Figure|2] illustrates that our wind model 
with a final mass-loss rate of 7xl0" 7 M o yr _1 is also able to 
reproduce the emission features due to amorphous silicates at 
9.7 and 18yum quite well. This successful model will be used 
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Fig. 3. Deconvolved 70 jim PACS image of R Aqr overplotted with 
contours from an ESO Multi-Mode Instrument (EMMI) 2007 observa- 
tion at 653 nm [Nn]. The contours clearly show a ring structure while 
the far-IR image renders only a part of this structure. 



in Sect. 14.1.31 to compare the dust temperatures resulting from 
the wind model to the corresponding ones estimated from PACS 
fluxes. 




Fig. 4. Map from the ratio of flux values between the deconvolved 
70 /mi and 160yt/m PACS images of R Aqr overplotted with contours 
from an EMMI observation at 653 nm [Nn]. The colour code repre- 
sents the flux ratio between 70 and 160//m. Both arms show gradients 
in the flux ratio, which is a signature of continuum dust emission. Fur- 
thermore, the highest flux ratio occurs at the ends of both arms. 



4.1 .3. Infrared arms in R Aqr: Jets, expanding nova, or ISM 
interaction? 

R Aqr was observed on the 16th of May 2011 with the PACS 
photometer of Herschel at 70/mi and 160/zm. Figure [T] shows 
the first successful attempts to image the far-IR emission around 
R Aqr. Both images reveal the same features. Since the 70/mi 
image offers a better resolution, however, it will be the base for 
the following discussion (unless otherwise noted). 

The far-IR circumstellar environment is basically charac- 
terised by two arms, one approximately in north-east direction 
and the other in south-west direction. In the 160yum band, the 
eastern arm suffers from the lower resolution and is only visible 
as a bright clump, whereas the western arm sustains its shape. 

The arms of R Aqr are exactly opposite at PA = 61° and 241° 
and seem to be connected with the binary system. However, ow- 
ing to the saturation of the central source, a look farther inside 
a region closer than « 18" from the star is not possible. The 
overall size of the IR emission structure around R Aqr is about 
95" x 35"0 20330 x 7490 AU). While both arms are aligned 
nearly symmetrically with respect to the central star, the light 
distribution is not. The western arm is brightest at its end, in the 
form of a light spot, whereas the eastern arm is clearly domi- 
nated by emission from the curved region. We now discuss three 
possible origins for the arms observed in the infrared. 



Bipolar jet? 

The first hypothesis naturally arises from the morphology of the 
IR emission and assumes th at the a rms trace the bipolar jet from 
the symbiotic system (Sect. 14.1.11) . This hypothesis is not well 
supported because of the discrepancy between the position an- 
gles of the IR arms (61° and 241°) and of the optical-UV-X jet 



(46° and 211°:lKellogg et al.l l2007l). Observations in the UV by 
Meier & Kafatos (1995) indicate that the direction of the bipo- 
lar wind is not stable due to the precession of the accretion disk. 
Although such a precession effect could perhaps account for the 
above discrepancy, we do not engage in discussing the possible 



association of the infrared arms with the jet, since the hypothesis 
cannot be validated with our observations. 



Nova nebula? 

The second hypothesis about the nature of the infrared arms thus 
stems from their good agreement with the structures seen for 
instance in [N n] 653 nm, as shown in Fig.[3j which correspond to 
th e equatorial plane of the outer hourglass structure first outlined 
bv lSolf & Ulri3hl (fl985h (see Sect. l4TTt Fl 

In the eastern arm, the optical emission is stronger on 
the northern half, which is moving away from the observer 
dHollis et al.ll 1999h . There it fits its far-IR counterpart well. In 
the southern part of the eastern arm, far-IR emission is missing, 
but the optical emission is also much weaker. Similar conclu- 
sions hold for the western arm: this time, it is the southern part 
of the ring that matches the IR emission (even the IR spot at 
the arm's end has a well-defined optical counterpart), whereas 
the northern part of the ring is missing IR emission. This close 
match between optical and far-IR features raises the question of 
the origin of the far-IR emission. There are two possibilities: ei- 
ther scattering of light from the central source by the gas or dust, 
or else continuum dust emission. In the first case, the good match 
between the optical and far-IR images is a natural consequence 
of the model. The high degree of polarizatio n observed from the 
nebula around R Aqr (lYudin & Evansll2002l) is in favour of this 
model. The second hypothesis requires the dust to be strongly 
coupled with the gas: it must either have formed in the nova 
ejecta traced by the gas seen in the optical image, or it must have 
been piled up by the ejecta in the Mira wind pre-existing to the 
explosion. 



2 The jets at PA = 46° and 211° ( Kellogg et alJ2007l) are clearly visible 
in the contour plot of Fig.[3]and do not match the far-IR emission. In the 
optical, the jets extend up to 20" north and 30" south from the star, well 
beyond the PACS far-IR emission, and have no IR counterpart, which is 
also true for the detached clump seen in the optical approximately 40" 
south of the star and for the (non-detached) clump 37" north. 
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To distinguish between the two models - light scattering or 
continuum dust emission - a map of the 70/mi to 160yum flux 
ratio was derived (Fig. [4]). Fluxes below 4<x of the background 
noise in either band were set equal to zero to prevent artefacts 
originating from background fluctuations. The part of the west- 
ern arm closest to the star suffers from artefacts in the 160yum 
band, which are not entirely eliminated by the 4cr clipping. The 
image clearly shows gradients in both arms, which are not ex- 
pected in the scattered-light scenario, since the flux ratio is a 
function of the scattering cross section only (provided that the 
scattering centres are the same throughout the arms) and is thus 
expected to be constant throughout the arms. The flux ratio im- 
age of Fig. |4] thus seems in favour of dust emission rather than 
light scattering. Even more surprising, the map reveals that the 
highest ratio occurs at the end of both arms, which are situated 
opposite each other much more accurately than in Fig. [3] At the 
end of the arms, the flux ratio is « 2.5 while at the northernmost 
part of the eastern arm, the ratio drops to « 1.1. To constrain 
the temperature from these flux ratios, it is necessary to specify 
the dust species in the arms. Assuming the dust to be constituted 
only of silicates with absorption coefficients Qabsjo = 3 x 10" 3 
and Qabs. 160 = 5.9 x 10" 4 at 70/mi and 160/zm, respectively 
(|Prainel ll985). the flux ratios observed in th e eastern arm corre- 
spond to temperatures of 37 K and 29 K (see lJorissen et alj201lL 
for details). Be cause astronom ical silicate has a very low opac- 
ity according to lDrainel(ll985l) . these values represent the lowest 
possible ones. 

The temperature of the eastern arm tip may now be compared 
with the predictions from the modelling at the same distance, 
since the model for photosphere and dusty outflow has been 
shown to describe the resulting SED of central source and inner 
wind region quite satisfactorily (Fig. 13 . The wind model yields a 
silicate dust temperature of 49 K at a distance of 9200 AU. These 
values compare reasonably well to the estimates based on the 
observational data (see above), and the successful application of 
the wind model to predict the dust temperature implies moreover 
that the dust is heated by the radiation from the central star rather 
than by some form of shock. 

The temperature gradient predicted by the wind model 
(Fig. [5} may now be used to de-project the geometry of the east- 
ern arm. Since the (logarithmic) gradient appears to stay con- 
stant for distances larger than 500 AU from the central star (i.e., 
well encompassing the location of the eastern arm), it is more 
robust than the temperature prediction itself (the dust tempera- 
ture of 37 K observed at a distance of 9200 AU differs somewhat 
from the 49 K value predicted by the model). The observed tem- 
perature ratio of 1.27 (= 37/29) along the eastern arm is thus 
equivalent to a (de-projected) distance variation by a factor of 
2.03 [= d(29 K)/d(37 K)]. Since the tip of the arm is its warmest 
part, it may be assumed that it is located closest to the star, by 
the factor 2.03 derived above, as compared to the coolest, north- 
ernmost point on the eastern arm. On the image, however, the 
corresponding (projected) distances are d pro j(29K) = 4300 AU, 
and J pr0 j(37 K) = 9200 AU, or a ratio of 0.467. This is in con- 
trast to the de-projected one, thus suggesting that the arm is 
strongly inclined onto the plane of the sky. To fix the ideas, we 
assume that the warmest tip of the eastern arm [at PA= 90°, cor- 
responding to dproj(37 (K)] lies in the plane of the sky (basically 
because the western and eastern tips are almost exactly oppo- 
site). With this assumption, we obtain <i pro j(37 K) = d(31 K) and 
d proj (29K)/d(29K) = cos/ = 0.467/2.03 = 0.23, or i = IT. 
Incidentally, this v alue is close to the 72 ° inclination of the or- 
bital plane found bv lSolf & Ulrichl (0~985). Therefore, the picture 
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Fig. 5. Tempe rature gradient delivered by the stationary wind model 
(see Sect. l4.1.2t . The dust temperatures of the eastern arm derived from 
the PACS band photometry (37 K and 29 K; blue triangles) correspond 
to distances of 21 500 AU and 43 700 AU in the wind model. 



emerging from our analysis for the arms is that of a non-circular 
structure located in the orbital plane. 

Interaction with the ISM? 

A third possible interpretation of the arm feature relates 
to its orientation with respect to the space motion of the 
star. Adop ting the proper motion from the SiO maser mea- 
surements (iKamohara et al.l 120101) and removing the solar 
motion (U, V, W ) <p = (11.10, 12.24,7.25) kms" 1 obtained by 
ISchonrich et al.l (l2010h . we found that the star moves with 
vlsr = 39.7 km s" 1 along PA = 145.1° with an inclination 
i = -40.2° with respect to the plane of the sky. In Fig. [T] 
the space motion is indicated by an arrow and its length gives 
the position of the star 2000 years in the future. The space 
motion thus appears to be almost perpendicular to the arms 
(PA = 61°). It is worth noting that a similar orthogonality has 
been found between the jet and the space motion of X Her 
(iJorissen et al.ll201 IT) . However, unlike the situation prevailing 
for X Her, the arm structure does not seem to be linked to the 
space motion, since ram pressure would bend the arms in the 
direction opposing the star motion, unlike what is observed for 
the eastern arm. Because of this, it is very unlikely that the ISM 
shapes the arm structure. 

As a conclusion to the discussion of these three possible 
models, the arms seem to be the IR counterpart of the expanding 
nebula associated with the outer hourglass observed in the 
optical and the emission is caused by radiation-heated dust 
grains. However, the temperature gradient within the eastern 
arm, along with the exactly opposing locations of the highest 
temperatures, seems to contrad ict the circul a r geo metry of the 
expanding ring proposed by ISolf & Ulrichl ([1985), since this 
model would imply if not a uniform temperature distribution, 
at least a mo r e arbi trary one than that observed. Although 
ISolf & Ulrichl (Il985l) clearly favoured the interpretation of an 
inclined ring ejecta, they also mentioned (their Sect. 4.3) the 
possibility that the nebula could be elliptical, based on kinematic 
arguments. This interpretation is supported by the temperature 
map of the Herschel observations. 
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Fig. 6. Evolutionary tracks from the Padova set (Ber telli et al.l l2008, 
2009) for 7=0.26 and Z=0.017, with masses between 1 and 8 Af , as in- 
dicated (in red) at the end of each track. The asymptotic giant branches 
are depicted with dashed (red) lines. The locati on of W Aql as de- 
rive d from the inte r ferom etric radii measured by iDanchi et al.l \ 19941) 
and Ivan Belle et al.1 dl997h are depicted by (red) triangles. 



4.2. WAql 

4.2.1. Background 

W Aql is a Mira v ariable star of spectral typ e S6.6 with a pe- 
riod of 490 days dFeast & Whitelockl I2000) at a distance of 
340 pc dGuan dalini & Busso ]j2008h . W Aql was found to be a 
binary star by Herbig (ll965l) . who noticed that the spectrum 
at minim um light acquir e d spec tral characteristics of an F5 or 
F8 star. iRamstedt et al.l (l201ll) recently presented an archive 
HST 435 nm image from 2004 that clearly shows the companion 
south-west of W Aql at a distance of 0'/47. This angular separa- 
tion corresponds to a projected distance of 160 AU at 340 pc. 

W Aql shows strong emission lines of SiO and HCN with 
a parabolic shape indicating that t he emission is optica lly thick 
(Bie ging et aDll998l) . In addition. IRamstedt et al.1 d2QQ9h found 
weak wings extending beyond the parabolic profile. They could 
be caused by recent mass-loss-rate modulations or asymmetric 
outflows. 

The mid-IR spectrum of W Aql is also dominated by emis- 
sion from amorphous silicate grains, with an additional contri- 
bution from amorphous alumina dust. Both emissions result in 
a peak position (for the blended silicate plus alumina bands) of 
about 10.5 /mi (lYang et al.ll2007l) . 

The circumstell ar e nvironment was stu died at 11.15yum by 
IDanchi et al] (1 19941) and lTatebe etaD(l2006l) . The former authors 
detected dust shells at 8 and 19 stellar radii while the latter found 
the dust to be unevenly distributed, with a clear excess eastward. 

4.2.2. Nature of the companion 

The suspicion of binarity ba sed on t he co mposite spectrum at 
minimum light discov ered by Herbig ( 1965]) was later confirmed 
by the HST images of lRamstedt et al.1 1 (1201 ll) . The log file of the 
available HST images is given in Table [5j Two sets of images 
separated by 1 1 yr are available and can be used to check for a 
possible relative motion of the two components. None has been 
detected, though, within the error bar set by the pixel size, which 
is 0Y04 for the WFPC, about two times larger than that of the 



Table 5. Images of W Aql retrieved from the HST archives, with col- 
umn labelled ID providing the HST program number. Column Am pro- 
vides the magnitude difference (m WAq i - m comp ) as estimated from aper- 
ture photometry. Columns labelled V and Vcomp provide the visual mag- 
nitudes of W Aql as derived from the AAVSO (for 1993) or ASAS (for 
2004) databases (Pojmanski 2002), and of its companion, respectively. 
The last column gives the measured separation 5. 



Instr. 
Filter 


Date 


ID 


Am 


V 


^comp 


6 
(") 


WFPC 
F487N 


1993-05-25 


3603 


+0.30 








WFPC 
F502N 


1993-05-25 


3603 


-0.99 


13.8 


14.8 


0.47 


ACS 
F435W 


2004-10-12 


10185 


-3.53 








ACS 
F606W 


2004-10-12 


10185 


-3.82 


11.1 


14.9 


0.47 



ACS image. At a distance of 340 pc, this pixel size corresponds 
to 13.6 AU. Assuming a circular orbit with a radius of 160AU 
(as derived from the angular separation of 0'.'47), the pixel size 
represents about 1.3 % of the orbital circumference. Over a 1 1 yr 
timespan, a motion would be detectable with the HST if the or- 
bital period were shorter than about 800 yr. 

A coincidental accord with a background source can also be 
excluded. The distance covered over 1 1 yr at a distance of 340 pc 
by an AGB star moving (rather slowly) with vlsr = 20 km s" 1 is 
more than a factor of 3 larger than the pixel size of the WFPC 
image (0V04 compared to 0V13 covered by the star). To de- 
rive the nature of the companion, aperture photometry has been 
performed on the two continuum images (obtained with filters 
F502N and F606W), and the results are listed in Table in the 
form of a magnitude difference with respect to W Aql. 

The magnitude of the Mira variable at the time of the HST 
observations was obtained from the AAVSO database (for the 
1993 observation) or from the ASAS database dPoj manskill2002l) 
(for the 2004 observation). It is remarkable that both epochs 
yield a consistent value for the visual magnitude of the com- 
panion (about 14.8), and that this magnitude is well below the 
magnitude of W Aql at minimum light (VwAqi,min ~ 14.0). This 
explains the absence of a flat minimum for the light curve, which 
would be caused by the companion becoming brighter than the 
variable star around minimum light. 

Therefore, adopting an apparent magnitude of 14.8 for the 
companion and a distance of 340 pc, the resulting absolute visual 
magnitude is My = 7.1, corresponding to a K4V main-sequence 
star with a mass of « 0.7 M Q This result co ntrasts with the F5 or 
F8 spectral type inferred by iHerbigl (Il965l) from the composite 
spectrum observed at minimum light; but the spectrum quality 
was quite poor. 

We now estimate the mass of W Aql to evaluate the or- 
bital period. IDanchi et all (11994 and Ivan Belle et al.l (1 19971) de- 
termined the angular diameter 6 and effective temperature of 
W Aql. Van Belle et al. (1997) observed the star at maxi- 
mum light and found 6 = 11.58 mas and = 2 230 K, which 
yields L = 18 590L© at a distance of 340 pc. IDanchi etaD 
(1994) observed the star at maximum and minimum light and 
found # max = 8.6 mas at r e ff 9max = 3000 K and # m i n = 7.2 mas 
at r e ff,min = 2500 K, corresponding to L max = 28 885L© and 
Lmi n = 9740 L Q . These values are widely spread and lead to a 
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Fig. 7. Deconvolved PACS images of W Aql at 70 /um (left) and at 160yt/m (right). The arrows indicate the direction of space motion and the 
position of the star 8000 yr in the future. 
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Fig. 8. Positional data of W Aql. The positions are taken from 9 
catalogues (see Tableland text) and cover the years 1913 to 2002. The 
(red) lines indicate the linear fit to the data points. The scatter is of the 
order of the angular separation of the binary system on the sky that is 
expected for a variability-induced movement. The positions of W Aql 
in the last century show a trend for a movement along PA = 38° with a 
space velocity v LS r = 21.6 + 4.1 kms -1 . 



stellar mass in the range 2 to 5 M© using the Padova stellar evo- 
lutionary tracks for a solar composition (Fig. |6). This yields a 
total system mass in the range 2.7 - 5.7 M . Adopting a circular 
orbit of radius 160 AU (0747 at 340 pc and i = 0), a minimum 
orbital period in the range 850 - 1230yr is found from Kepler's 
third law. This estimate agrees with that based on the absence of 
any relative motion on the HST image and has to be seen as the 
lowest bound, since the de-projected radius of the orbit can be 
much larger than 160 AU. 



4.2.3. Space motion 

Unfortunately, the space motion of W Aql is not clearly defined, 
as shown in Table [3j which collects proper-motion values from 
various catalogues. Unlike the situation prevailing for the proper 
motions of R Aqr, also listed in Table [3] the proper motion of 



W Aql varies significantly from one catalogue to the other, with 
no trend with the mean epoch spanning almost a century. 

This behaviour very likely results from the combination of 
the Mira pulsation and the presence of a companion. During 
the variability cycle of the Mira, the photocenter of the binary 
system moves along the path that joins the two stars (separated 
by 0.47"), according to the respective contribution of the two 
components to the total light (according to AAVSO, the Mira 
variable varies from V = 7 to 14, while the companion has V ~ 
14.8 as shown above). The proper motion, of the order of a few 
tens of milliarcseconds per year, can therefore not be reliably 
derived, since the individual position measurements on which it 
relies are varying over a range of about 400 milliarcseconds , with 
a periodicity of 490 d. This behaviour is denoted as "variability- 
induced mover" (VIM) in the c ontext of the Hipparcos catalogue 
(see also lPourbaix et al.ll2003l) . 

Nevertheless, on a century-long time span, one may hope 
that the long-term proper motion will dominate over this short- 
term variability-induced motion. For this reason one can try to 
derive the proper motion from the positional data of the seven 
catalogues in Table[3](see references within) for which an obser- 
vation epoch is stated. The SAO Star Catalog has to be excluded 
because W Aql's position is far off from that of the other cat- 
alogues and seems erroneous. Instead, the positions from th e 
HST Guide Star Catalog (obs. epoch 1987 jLasker et al.lll996l) . 
the Tycho-2 Catalog (obs. epoch 1991 :lH0get a l. 2000) and 
the USNO-B Catalog (obs. epoch 2000: lMonet et al.ll2003h are 
added. They are not listed in Table [3]because they do not feature 
proper motion values for W Aql. The result are nine data points 
that cover nearly one century (1913-2002). Figure [8] shows the 
right ascension and declination values with respect to the obser- 
vation epoch. The diagrams clearly show a trend for the stellar 
motion of W Aql. Linear curve fitting was performed to ob- 
tain the proper motion from the slopes of the curves. The re- 
sulting proper motion values are fi* a = 11.8 ± 4.1 mas yr -1 and 
[is = 4.5 ± 2.1 mas yr -1 for right ascension and declination, re- 
spectively. These values are different from all stated in Table [3] 
but reflect observations of the stellar position over a timespan of 
nearly one century, essentially longer than the observations for 
every catalogue in the table. In addition, the direction of motion 
is largely supported by the IR observations of the CSE, which 
are discussed in Section 14.2.41 

Although the scatter around the linear fit seems large at first 
sight, it is « 2 x 10" 4 degrees or « 077, which is about the 
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Fig. 9. Left panel: Contour plot of the deconvolved 70 /dm PACS image of W Aql. Two Archimedean spirals matching the outline of the IR 
emission are recognizable. The (black) dotted spiral has one winding to cover the emission with a separation distance 6 = 52" while the (white) 
dashed spiral has two windings with a separation distance 6 = 26". Right panel: the same image overplotted with an ellipse fitted to the IR 
emission. The axes lengths of the ellipse are 69" and 59", respectively, with the major axis oriented along PA= 35°. 



separation of the binary system (0'/5). This is exactly what is 
expected for a variability-induced movement, whose amplitude 
must be of the order of the angular separation on the sky. 

Adopting a distance of 340 pc ( Guandali ni & Bussol 
2008 j), a radial velocity V r = -20.80 ± 3.10kms" 1 
dKharchenko et al.1 [2007b a nd the solar mo ti on (U V, W) Q = 
(11.10, 12.24, 7.25)kms" 1 dSchonrich et al.1 l2010h . the re- 
sulting proper motion values are ji* a = 7.9 ± 4.1 mas yr" 1 
and fi$ = 10.0 + 2.1 mas yr" 1 for right ascension and dec- 
lination, respectively. The space velocity of W Aql is 
vlsr = 21.6 ± 4.1 km s" 1 at PA= 38° with an inclination 
i = -17° with respect to the plane of the sky. However, due to 
the large uncertainties of the positions, the results thus obtained 
may include large errors. 

4.2.4. Binary vs. ISM interaction 

W Aql was observed with PACS on October 2011. Figure □ 
presents the 70//m and 160/mi images. Since they exhibit the 
same features, the discussion below relies on the 70 jim image. 
In the IR, the CSE is at first sight a ty pical fermata case and was 
classified as such by ICox et all (120121) . The elliptical envelope 
shape is to be expected from the interaction of the stellar wind 
with the oncoming ISM and the orientation of the envelope with 
its symmetry axis at PA « 35° (with the star offset by 8" though) 
nicely fits the space motion of the star along PA = 38° (see right 
panel of FigO. 

A feature that does not directly fit with this interpretation 
though is the brighter region east of the star, at a distance of 
42". This feature suggests a bow shock but its location and the 
space motion/symmetry axis of the ellipse do not match: the 
bright region eastward extends from PA = 55° to 140° while 
the stellar motion and the CSE have their symmetry axis along 
PA « 35°. However, it is imaginable that a strong ISM flow 
blowing from the east modifies the CSE morphology and creates 
a fermata. Another possible scenario is an ISM flow from the 
south that pushes material to the east. A feature supporting 
this is the bending of the IR emission « 28" south of the star 
towards the star. Such ISM flows combine d with the prope r 
motion of the star were proposed for a Ori by lUeta et al.l (120081) . 
To infer the direction of the ISM flow, it is necessary to know 
the three-dimensional orientation of the bow shock, which is 



virtually impossible since the inclination of the bow shock 
cone has almost no effect on the (obs erved) projected outline 
dMac Low et al.lll99ll: ICox et alJ l2012h . Thus, the hypothesis 
of an ISM flow that interacts with the CSE at the bright region 
could be adequate, but it is hardly verifiable. 

It is hard to imagine, however, that an ISM flow that cre- 
ates such a strong shocked interface does not modify the overall 
symmetry of the CSE, which is aligned with the (quite slow) 
space motion of the star. The binary nature of W Aql made 
us check whether a companion-induced Archimedean spiral 
could fit the outline of t he CSE, as s een around RW LMi 
(Mauron & Huggins 2006) and LL Peg dDinh-V.-Trung & Liml 
2009). Mastrodemos & Morris (1999) showed in their hydro- 
dynamic simulations that for large binary separations, such 
spirals even ap pear as broken a rcs or concentric shells when 
seen edge-on (iMayer et~ai1l201 ll) . The fact that such structures 
are not restricted to face-on orbits increases the probability to 
observe them. 

Given the orbital period P > lOOOyr and the wind velocity 
v w = 20kms _1 (Section 14.2.21 and Table [TT). an arm separation 

p = v w xP> 12'/5 (5) 

is estimated. Within this range, we indeed find two solutions, as 
drawn on the left panel of Fig. [51 with arm separations of 26" 
or 52", the former being about twice the estimated (minimum) 
winding separation. On the western side, the shell around W Aql 
exhibits three equidistant sharp ridges that are quite clearly de- 
lineated by the 26" Archimedean spiral. That spiral is broken on 
its southern side, however. An orbital period of 2000 yr is thus 
required to match the Archimedean spiral, which is certainly not 
impossible, given the uncertainties associated with the distance 
and the orbital inclination. 

In addition, as indicated bv lRaga et all (l201lh . the starting 
point of the spiral structure ought not to be at the centre of the 
binary, but is offset by a distance 



where r or b is the orbital separation, and v or b is the orbital veloc- 
ity. Within this distance, the wind propagates shock-less. For 
the W Aql binary, the estimated orbital separation is 160AU 
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and the orbital velocity 2.4 km s" 1 (assuming an orbital period 
of 2000 yr). This gives a shock- free radius of 1330 AU or 4". 

In summary, there are good arguments to support both hy- 
potheses of a binary-induced spiral and of a wind-ISM inter- 
action. Although the Archimedean spiral nicely fits the CSE, 
it cannot explain the bright region in the east. Similarly, an 
ISM flow alone cannot explain the CSE shape that follows an 
Archimedean spiral. Therefore, it is suggested that an ISM flow, 
preferentially from east or south direction, interacts with a part 
of the spiral and forms the bright region. A similar combination 
of orb ital spiral and ISM interaction was already reported around 
o Cet dMaver et al.ll201lh . 

5. Conclusions 

We have analysed the properties of the 18 targets among the 78 
AGB stars and RSGs in the M ESS samp l e tha t are suspected 
of being binaries. Based on the Co x et al. classification 
of the circumstellar morphology into five groups, the binary ob- 
jects tend to be dominantly of irregular type (no imprint of the 
wind-ISM interaction on the circumstellar morphology). In con- 
trast, the ring-type category contains only one binary star, which 
suggests that the interaction of the companion with the stellar 
wind decisively alters the spherical symmetry of the wind. How- 
ever, despite the small number of binary stars in the sample, 
we showed that the null hypothesis of a similar binary fraction 
among the irregular shells and among the total sample can be 
rejected at a 98% confidence level. 

Among the 17 binary AGB stars of the sample, a detailed 
analysis of the PACS images at 70/mi and 160 /mi was per- 
formed for the two objects R Aqr and W Aql. 

R Aqr, a symbiotic star, shows two arms extending over a 
region 95" x 35" wide. The arms are exactly opposite at PA 
= 61° and 241°, b ut their position doe s not fit the X-ray jets 
(PA 41° and 211°; iKellogg et al]l2007l) . However, a contour 
plot drawn from an ESO Multi Mode Instrument (EMMI) op- 
tical image shows that the IR ar ms partly fit the ring nebula de- 
scribed bv lSolf &Ulrichl dl985). A map of the flux ratios be- 
tween 70yum and 160yum, which is equivalent to a temperature 
map, reveals a gradient in both arms with the highest tempera- 
ture at their ends. This result dismisses a circular-ring structure, 
which would not show any temperature gradient in the arms, and 
favours an elliptically- shaped nebula. In addition, the SED of 
R Aqr (constructed with observational data compiled from the 
literature) was fitted with the help of a combined model for pho- 
tosphere and dusty outflow, the latter having a mass-loss rate of 
M^7xl0~ 7 M o yr _1 . A silicate dust temperature of 37— 49 K 
for the tip of the eastern arm could be estimated from the PACS 
photometry as well as from the modelling of the stellar wind. 
The temperature gradient delivered by the wind model results 
in an inclination of the arms by 77°, which agrees well with the 
inclination of the orbital plane (72°: ISolf & Ulrichl 19851) . There- 
fore, the arms are part of a non-circular structure located in the 
orbital plane. 

For the second object, W Aql, Ramsted t et al.l (1201 ll) pre- 
sented an HST 435 nm-image of the companion, which was pre- 
viously detected spectroscopically as an F5 or F8 star (Herbig 
1 19651) . The quality of Herbig 's spectrum is quite poor, how- 
ever. Aperture photometry on archive HST images at other 
wavelengths (502 and 606 nm) and light-curve data reveal that 
the companion has an absolute magnitude My = 7.1, corre- 
sponding to a K4 main-sequence star with a mass of « 0.7 M . 
With a mass for W Aql between 2 and 5 M Q (as determined 
from HRD track matching on luminosity and angular diameter 



data bv lDanchi etaD (1 19941) and lvan Belle et al] (Il997h ). one de- 
rives an orbital period of at least lOOOyr assuming a separation 
of 160 AU. The outline of W Aql's far-IR emission can be re- 
produced by an Archimedean spiral that has either one or two 
windings, with an arm separation of 52" or 26", respectively. A 
separation of 26" would imply that the orbital period is about 
twice as large as the inferred lOOOyr. However, this value for 
the period assumed an orbit seen face-on, but the binary sepa- 
ration (and hence the orbital period) may be much larger if the 
orbital plane is inclined with respect to the plane of the sky. The 
two-arm Archimedean spiral thus appears as a viable solution. 
But the CSE also shows some features, especially a bright re- 
gion east of the star that can only be explained by the interaction 
with the ISM. Although the symmetry axis of the IR emission is 
along the space motion at PA= 36° and differs from the bright 
region at PA 55-140°, an ISM flow from the east or south-east 
directions is a valid explanation for this feature. Thus the wind 
of W Aql is most probably shaped both by the companion, which 
produces the spiral, and by the ISM, which is responsible for the 
spiral brightening upstream of the ISM flow. In a forthcoming 
paper, the binaries 6 Aps and n l Gru will be analysed along sim- 
ilar guidelines. 
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